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ABSTRACT 

We observed near-Earth asteroid (162421) 2000 ET70 using the Arecibo and Goldstone radar sys- 
tems over a period of 12 days during its close approach to the Earth in February 2012. We obtained 
continuous wave spectra and range-Doppler images with range resolutions as fine as 15 m. Inversion 
of the radar images yields a detailed shape model with an effective resolution of 100 m. The asteroid 
has overall dimensions of 2.6 km x 2.2 km x 2.1 km and a surface rich with kilometer-scale ridges 
and concavities. It is a principal axis rotator and spins in a retrograde manner with a spin period of 
8.96 ± 0.02 hours. In terms of gravitational slopes evaluated at scales of 100 m, the surface seems 
mostly relaxed with over 99% of the surface having slopes less than 30° , but there are some outcrops 
at the north pole that may have steeper slopes. Our precise measurements of the range and veloc- 
ity of the asteroid, combined with optical astrometry, enables reliable trajectory predictions for this 
potentially hazardous asteroid in the interval 460-2813. 

Subject headings: Near-Earth objects; Asteroids; Radar observations; Orbit determination 



1. INTRODUCTION 

Radar astronomy is arguably the most powerful Earth- 
based technique for characterizing the physical properties 
of near-Earth asteroids (NEAs) . Radar observations rou- 
tinely provide images with decameter spatial resolution. 
These images can be used to obtain accurate astrome- 
try, model shapes, measure near-surface radar scattering 
properties, and investigate many other physical proper- 
ties (e.g. sizes, spin states, masses, densities). Radar ob- 
servations have led to the discovery o f asteroids exhibit- 
ing non-principal axis rotation (e.g., Ostro et al."1995^ 
Benncr ct al. 2002), binary an d triple NEAs (e.g ., Margot 
et al. 2002; Ostro ct ah 2006; Sh epard et al.||2006j INoiaiTI 
et al. _2008j Brozovic et al. 2011J , and c ontact binary as- 
teroids ( ag.t Hudson and Ostrojl 994 Bc nner et al.|2006" 



Brozovic et aT 2010[ ). Radar-derived shapes and spins 



have been used to mvestigate various physical processes 
(Yarkovsky, YORP, BYORP, tides, librations, e tc.) that 
are important to the ev olution of NEAs (e.g., [Ches le^j 



et aH2003| [Nugent et al 2012; L owry et al.l20 07: 'I'aylor 
et al.l|2007| [Ostro et al.||2006, .Scheeres et al.||2006, .Tay- 



lor and Margot|2011[|Fang et al.|2011||Fang and Margot 
2012| ; 

Here we present the radar observations and detailed 
physical characterization of NEA (162421) 2000 ET70. 
This Aten asteroid (a=0.947 AU, e=0.124, i=22.3°) was 
discovered on March 8, 2000 by the Lincoln Near-Earth 
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Asteroid Research (LINEAR) program in Socorro, New 
Mexi co. Its a.bsolut e magnitude was reported to be 
18.2 (Whiteley 2001) which, taken at face value, sug- 
gests a diameter between 0.5 and 1.5 km (assuming an 
optica l albedo between 0.4 and 0.04). [Alvarez et al. 
(20121 obtained a hghtcurve of the asteroid during its 
close approach to the Earth in February 2012. They re- 
ported a lightcurve period of 8.947 ± .001 hours and a 
lightcurve amplitude of 0.60 ± 0.07. [Whiteleyj (|2001 1 
classified 2000 ET70 as an X-type asteroid whicn is a 
degener ate group of asteroids consisting of E, M, and P 
classes ( Tholen|1984 ). Mike Hicks (personal communica- 
tion) indicated that his observations were best matched 
by a C-type or possibly E-type. One of us (EH) cl assified 
it as Xk in the taxonomic system of Bus-DeMeo ( [DeMeo 
eral][2009t . 



2. OBSERVATIONS AND DATA PROCESSING 

We observed 2000 ET70 from February 12, 2012 to 
February 17, 2012 using the Arecibo S-band (2380 MHz, 
13 cm) radar and from February 15, 2012 to February 
23, 2012 using the Goldstone X-band (8560 MHz, 3.5 
cm) radar. The asteroid moved ~74 degrees across the 
sky during this time and it came closest to Earth on 
February 19 at a distance of ~0.045 Astronomical Units 
(AU) (Fig.jl]). We observed it again in August 2012 when 
it made another close approach to Earth at a distance of 
-0.15 AU. 

Radar observing involved transmitting a radio wave for 
approximately the round-trip light-time (RTT) to the as- 
teroid, ~45 seconds at closest approach, and then receiv- 
ing the echo reflected back from the target for a compa- 
rable duration. Each transmit-receive cycle is called a 
run. On each day we carried out runs with a monochro- 
matic continuous wave (CW) to obtain Doppler spec- 
tra, followed by runs with a modulated carrier to obtain 
range-Doppler images. Table [l] summarizes the CW and 
range-Doppler imaging runs. 
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TABLE 1 

Radar Observations of (162421) 2000 ET70 



Tel 


UT Date 


MJD 


Eph 


RTT 


PTX 


Sr 


5f 


Code 


Start-Stop 


Runs 




yyyy-mm-dd 








(kW) 


(m) 


(Hz) 




hhmmss-hhmmss 




A 


2012-02-12 


55969 


s41 


69 


828 


cw 


0.167 


none 


082751-083755 


5 








s41 






15 


0.075 


65535 


084247-102947 


48 








s41 






15 


0.075 


8191 


105318-110955 


8 


A 




Oo J I U 




D4 


OOU 


cw 




none 


UoiiUD-UoZoUo 











s43 






15 


0.075 


65535 


083034-105326 


54 


A 


2012-02-14 


55971 


s43 


60 


811 


cw 


0.196 


none 


075956-080443 


3 








s43 






15 


0.075 


65535 


080640-101945 


59 


A 


2012-02-15 


55972 


s43 


56 


785 


cw 


0.213 


none 


074302-074729 


3 








s43 






15 


0.075 


65535 


080301-080918 


4 








s47 






15 


0.075 


65535 


081638-093235 


42 








s47 






15 


0.075 


65535 


093749-100946 


18 


G 


2012-02-15 


55972 


s43 


56 


420 


cw 




none 


085515-090327 


5 








s43 






75 


1.532 


255 


091748-093320 


9 








s45 






37.5 




255 


094624-115924 


73 








s45 






37.5 




255 


094624-115924 


5 


A 


2012-02-16 


55973 


s51 


52 


760 


cw 


0.227 


none 


073418-073830 


3 








s49 






15 


0.075 


65535 


074056-074652 


4 














cw 




none 


074838-075106 


2 














10 


U.U ( o 


DOOoO 


U ( Oozo-uyoolO 


DZ 


G 


2012-02-16 


55973 


s49 


52 


420 


cw 




none 


091546-092330 


5 








s49 






75 


1.532 


512 


095639-100239 


4 








s49 






37.5 


0.488 


none 


113521-114633 


7* 








s4y 






1 

io 


l.U 


none 


izi004-iz4 1 lU 


io 








s49 






15 


1.0 


none 


131001-132809 


11* 








s49 






15 


1.0 


none 


132904-152931 


70* 


A 


2012-02-17 


55974 


s53 


49 


775 


cw 


0.244 


none 


073800-074157 


3 








SOO 






1 ^ 
10 


U.U ( o 


DOOoO 


U ( 44ioD-Uo4ooy 


4U 


G 


2012-02-17 


55974 


s53 


49 


420 


cw 




none 


070553-071310 


5 








s53 






75 


1.532 


512 


074255-080001 


11 








SDo 






O 1 .O 


U.y ( i 


none 


Uoiuo i-izz4iy 


-| CO* 


G 


2012-02-18 


55975 


s53 


47 


420 


CW 




none 


070551-071250 


5 








s53 






75 


1.532 


512 


073604-075052 


10 








s53 






37.5 


0.977 


none 


080115-083144 


20* 








SDo 






O i .O 


U.y ( i 


none 


nCific;-! HQ/I c;0/1 


D 


G 


2012-02-19 


55976 


s55 


46 


420 


CW 




none 


070551-071241 


5 








s55 






75 


1.532 


512 


072155-973625 


10 








SOO 






O / .0 


U.y 1 1 


none 


U 1 4DiiS-i i ioOo 


ioD 








s55 






37.5 


0.977 


none 


114410-130706 


55* 


G 


2012-02-20 


55977 


s57 


45 


420 


cw 




none 


071051-071742 


5 








so / 






ol .0 


U.y 1 1 


none 


in'51/IO 11 C^O 

iUoi4y-i ioooy 


/1 1 * 

41 








s57 






37.5 


0.977 


none 


121608-121954 


3* 


G 


2012-02-22 


55979 


s59 


48 


420 


cw 




none 


094825-095542 


5 








s59 






75 


0.957 


512 


100253-102135 


12 








s59 






37.5 


0.977 


none 


103251-103514 


2* 








s59 






37.5 


0.977 


none 


103606-104502 


6* 


G 


2012-02-23 


55980 


s59 


50 


420 


cw 




none 


083310-084053 


5 








s59 






75 


0.957 


512 


084948-090941 


12 








s59 






37.5 


0.977 


none 


092046-105520 


55* 


A 


2012-08-24 


56163 


s72 


153 


721 


cw 




none 


151051-151825 


2 














cw 




none 


154651-163117 


9 














2 


2 


8191 


163659-174101 


13 


A 


2012-08-26 


56165 


s72 


157 


722 


cw 




none 


150424-162038 


15 



Note. — The first column indicates the telescope: Arecibo (A) or Goldstone (G). MJD is the modified 
Julian date of the observation. Eph is the ephemeris solution number used (Section [sj. RTT is the 
round-trip light-time to the target. PTX is the transmitter power. 5r and Sf are the range and Doppler 
resolutions, respectively, of the processed data. Code is the length of the pseudo-random code used. 
The timespan of the received data are listed by their UT start and stop times. Runs is the number of 
transmit-receive cycles during the timespan. An asterisk (*) indicates chirp runs. 
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Fig. 1.— Distance of 2000 ET70 from Earth (top) and sub-radar 
point (SRP) latitude (bottom) during the primary observing pe- 
riod. SoUd and dashed Hnes in the bottom panel show SRP lat- 
itude for the prograde and retrograde spin vectors, respectively 
(Section [6]l. 

For CW runs a carrier wave at a fixed frequency was 
transmitted for the RTT to the asteroid. The received 
echo from the asteroid was demodulated, sampled, and 
recorded. A fast Fourier transform (FFT) was applied 
to each echo timeseries to obtain the CW spectra. The 
total frequency extent (A/) or bandwidth (BW) of the 
CW spectra is equal to the sampling frequency (fg) or 
the reciprocal of the sampling period {Pg)- 



A/ = /. = 



1 



The spectral resolution (Sf) is given by: 



Sf 



Af 

n 



(1) 



(2) 



where n is the FFT length. Finer spectral resolution 
can be achieved by increasing n, however the signal- 
to-noise ratio (SNR) in each frequency bin decreases 
as Xj^pn. The finest possible resolution that can be 
achieved is limited by the number of samples obtained 
in one run. If we were recording for the full duration 
of the RTT, the number of samples obtained in one run 
would be n = RTT x /g, and the resolution would be 
(5/ — 1/RTT. In reality we cannot transmit for a full 
RTT as it takes several seconds to switch between trans- 
mission and reception, so the finest possible resolution 



), where ^switch is the switching 



is (5/ = 1/(RTT - i^witch 
time. 

For range-Doppler imaging two different pulse com- 
pression techniques were used to achieve fine range reso- 



lution while maintaining adequate SNR (Peebles 2007). 
Pulse compression is a signal conditioning and processmg 
technique used in radar systems to achieve a high range 
resolution without severely compromising the ability to 
detect or image the target. The range resolution achiev- 
able by a radar is proportional to the pulse duration, or, 
equivalently, inversely proportional to the effective band- 
width of the transmitted signal. However, decreasing the 
pulse duration reduces the total transmitted energy per 
pulse and hence it negatively affects the ability to detect 
the radar target. Pulse compression techniques allow for 
the transmission of a long pulse while still achieving the 
resolution of a short pulse. 

2.1. Pulse compression using binary phase coding 
(BPC) 

We used binary codes to produce range-Doppler images 
of the asteroid at both Arecibo and Goldstone. In this 
technique we modulated the transmitted carrier with a 
repeating pseudo-random code using binary phase shift 
keying (BPSK). The duration of the code is called the 
pulse repetition period (PRP). The effective bandwidth 
of the transmitted signal (-Boff ) is given by l/Tbaud, where 
Tbaud is the duration of each symbol of the code. For 
each run we transmitted for approximately the RTT to 
the asteroid, followed by reception for a similar duration. 
The received signal was demodulated and then decoded 
by cross-correlating it with a replica of the transmitted 
code. The images span a range (Ar) given by: 

Ar = ^PRP, (3) 

and their range resolution {5r) is given by: 

(Jr = ^Tbaud = (4) 

where c is the speed of light. Each baud within the code 
eventually maps into to a particular range bin in the 
image. Resolution in the frequency or Doppler dimension 
was obtained in each range bin by performing a FFT on 
the sequence of returns corresponding to that bin. The 
total frequency extent (A/) or bandwidth (BW) of the 
image is equal to the pulse repetition frequency (PRF) 
or the reciprocal of the PRP: 



A/ = PRF 



PRP 



(5) 



The frequency resolution {5f) of the image depends on 
the FFT length (n) as follows: 



5f = 



A/ 



(6) 



The RTT dictates the finest frequency resolution achiev- 
able, similar to the situation with CW spectra. 

2.2. Pulse compression using linear frequency 
modulation ( Chirp ) 



We used a linear frequency modulation technique ( Pee- 
bles 2007; Margot 2001) to produce range-Doppler im- 
ages of the asteroid at Coldstone. Chirp waveforms allow 
us to maximize the bandwidth of the transmitted signal 
and to obtain better range resolution than that available 
with BPC waveforms. They are also less susceptible to 
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degradation due to the Doppler spread of the targets. 
Finally, they are more amenable to the application of 
windowing functions that can be used to trade between 
range resolution and range sidelobe level (Margot 2001). 

In this technique the carrier was frequency modulated 
with a linear ramp signal. The resultant signal had a 
frequency that varied linearly with time from loq — uj to 
ujQ+uj, where ojq is the carrier frequency and BeS — 2uj is 
the effective bandwidth of the signal. The resultant sig- 
nal is called a chirp. A repeating chirp was transmitted 
with 100% duty cycle for the duration of the round-trip 
light-time to the asteroid. The time interval between the 
transmission of two consecutive chirps is the PRP. For 
our observations we used PRPs of 125 fis and 50 iis for 
chirps with BcS=2 MHz and Bcff=5 MHz, respectively. 
The received signal was demodulated and range com- 
pression was achieved by cross-correlating the echo with 
a replica of the transmitted signal. The range extent 
and the range resolution of the chirp images are given 
by Eq. ([s]) and Eq. Q, respectively. Resolution in the 
frequency or Doppler dimension was obtained as in the 
binary coding technique. The bandwidth of the image is 
given by Eq. ^ and the frequency resolution is given by 
Eq. ^. 

3. ASTROMETRY AND ORBIT 

A radar astrometric measurement consists of a range or 
Doppler estimate of hypothetical echoes from the center 
of mass (COM) of the object at a specified coordinated 
universal time (UTC). In practice we used our measure- 
ments of the position of the leading edge of the echoes 
and our estimates of the object's size to report prelim- 
inary COM range estimates and uncertainties. We re- 
fined those estimates after we obtained a detailed shape 
model (Section [?]). We measured Doppler astrometry 
from the CW spectra. We reported 9 range estimates 
and 1 Doppler estimate during the course of the observ- 
ing run and computed the heliocentric orbit using the 
JPL on-site orbit determination software (OSOD). The 
ephemeris solution was updated each time new astromet- 
ric measurements were incorporated (Table [T]). 

At the end of our February observing campaign, we 
were using ephemeris solution 59. A final range mea- 
surement, obtained during the asteroid's close approach 
in August 2012, was incorporated to generate orbit so- 
lution 74. After the shape model was finalized, we up- 
dated the orbit to solution 76 by replacing the prelim- 
inary February astrometric measurements with shape- 
based astrometric measurements. Table [2] lists the final 
radar astrometric measurements. 

Table [s] lists the best fit orbital parameters (solution 
76) generated using 18 range measurements and 316 op- 
tical measurements. The optical measurements span 
February 1977 to December 2012. However we assigned 
20 arcsecond uncertainties to the two precovery observa- 
tions from the 1977 La Silla-DSS plates, effectively re- 
moving their contribution to the fit and reducing the 
optical arc to the interval 2000-2012. The 1977 obser- 
vations were from a single, hour-long, trailed exposure, 
and appear to have been reported with a ~45 s tim- 
ing error, consistent with the measurers' cautionary note 
"Start time and exposure length are uncertain" (MPEC 
2000-L19). 

The orbit computation is reliable over a period from 



the year 460 to 2813. Beyond this interval either the 3(t 
uncertainty of the Earth close approach time exceeds 10 
days or the Sti uncertainty of the Earth close approach 
distance exceeds 0.1 AU. Here encounter distances of less 
then 0.1 AU are considered as close approaches. The 
current Minimum Orbit Intersection Distance (MOID) 
with respect to Earth is 0.03154 AU, making 2000 ET70 
a potentially hazardous asteroid (PHA). 

4. RADAR SCATTERING PROPERTIES 

We transmitted circularly polarized waves and used 
two separate channels to receive echoes having the same 
circular (SC) and the opposite c ircular (OC) polarization 
as that of the transmitted wave ( Ostro|1993 ). Reflections 
from a plane surface reverse the polarization of the inci- 
dent waves and most of the echo power is expected in the 
OC polarization. Echo power in the SC polarization is 
due to multiple reflections or reflections from structures 
with wavelength-scale roughness at the surface or sub- 
surface. A higher ratio of SC to OC power therefore in- 
dicates a greater degree of near-surface wavelength-scale 
roughness or multiple scattering. This circular polariza- 
tion ratio is often denoted by fic- We measured ^iq for 
all the Arecibo spectra shown in Fig. [2] and computed an 
average value of fj,c = 0.21 ± 0.02, where the error bar is 
the standard deviation of the individual estimates. Ob- 
served ratios for individual spectra deviate no more than 
0.03 from the average. This ratio is lower than that for 
the ma jority of NEAs with known circular polarization 
ratios (Benner et"al]|2008[ ) suggesting that 2000 ET70 
has a lower than average near-surface roughness at 10 
cm scales. 

In our modeling of the shape of the asteroid (Sections|6] 
and[7|, we used a cosine law to represent the radar scat- 
tering properties of 2000 ET70: 



da 
dA 



i?(C + 1) (cos a) 



2C 



(7) 



Here a is the radar cross section, A is the target sur- 
face area, R is the Fresnel reflectivity, C is a parame- 
ter describing the wavelength-scale roughness, and a is 
the incidence angle. C is related to the wavelength-scale 
adirectional root-mean-square (RMS) slope Sq and angle 
of the surface by So = tan (^j-ms) = C~^/^ . Larger 



values of C represent more specular scattering (Mitchell 
eralin[996). 



5. RANGE AND DOPPLER EXTENTS 

The range extent of the object in the radar images 
varies between ~600 m and ^-^1700 m (Fig. 8]), suggest- 
ing that the asteroid is significantly elongateoT^In most of 
the images two distinct ridges that surround a concavity 
are clearly visible (e.g.. Fig. [8| images 8-12 and images 
30-37). In images where the ridges are aligned with the 
Doppler axis, they span almost the entire bandwidth ex- 
tent of the asteroid (e.g., Fig.Jsj images 11 and 34). If 
the concavity is a crater then these ridges could mark its 
rim. At particular viewing geometries, the trailing end 
of the asteroid exhibits a large outcrop with a range ex- 
tent of ^250 m (e.g.. Fig. |8| images 11). These features 
suggest that the overall siirface of the asteroid is highly 
irregular at scales of hundreds of meters. 

For a spherical object, the bandwidth {B) of the radar 
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Fig. 2. — Arecibo CW spectra of 2000 ET70 arranged in chronological order from left to right and top to bottom. They were obtained 
on February 12-17 (MJDs 55969.35, 55970.34, 55971.33, 55972.32, 55973.32, 55974.32, respectively). The frequency resolution of each 
spectrum is given at the top right of the corresponding panel. Each spectrum was produced by choosing a frequency resolution that allowed 
for the incoherent sum of 10 independent spectra, or looks, per run, and by summing over 3 runs, resulting in a total number of 30 looks. 
The solid and dashed lines are observed OC and SC spectra, respectively. The dot-dashed line shows the corresponding synthetic OC 
spectra generated using our shape model (Section[7]l. The circular polarization ratio (/^c) values in chronol ogical order are .18. 0.20, 0.24, 
0.21, 0.21, and 0.24 all of which have uncertainties of 5%, where the uncertainty is computed according to |Ostro et al.| l |l983| |. 



TABLE 2 

Radar astrometric measurements of 2000 ET70 



Date (UTC) 


Range 1-a 


Uncertainty 


Observatory 


yyyy-mm-dd hh:mm:sE 


i ixs 


IXS 




2012-02-12 08:50:00 


67220894.76 


0.5 


A 


2012-02-12 10:03:00 


66974197.50 


0.5 


A 


2012-02-13 09:11:00 


62452848.06 


0.5 


A 


2012-02-13 10:00:00 


62298918.42 


0.5 


A 


2012-02-14 09:28:00 


58061903.74 


0.5 


A 


2012-02-14 10:06:00 


57953990.18 


0.5 


A 


2012-02-14 10:16:00 


57925824.25 


0.5 


A 


2012-02-15 08:30:00 


54319877.92 


0.5 


A 


2012-02-15 09:20:00 


54206358.81 


3.0 


G 


2012-02-15 09:47:00 


54126135.16 


0.5 


A 


2012-02-16 08:23:00 


50979848.61 


0.5 


A 


2012-02-16 09:30:00 


50840187.02 


0.5 


A 


2012-02-17 08:05:00 


48327621.88 


0.5 


A 


2012-02-17 08:43:00 


48267222.55 


0.5 


A 


2012-02-18 07:40:00 


46481930.28 


2.0 


G 


2012-02-19 07:30:00 


45481292.56 


2.0 


G 


2012-02-20 10:10:00 


45483470.81 


2.0 


G 


2012-08-24 17:08:00 


153139162.28 


3.0 


A 


Note. — The first 


column indicates the 


coordinated 


universal time 



(UTC) of the measurement epoch. The second column gives the range 
measurements as the RTT to the asteroid in microseconds (/is). The third 
column lists the la range uncertainty. The fourth column indicates the 
radar used to make the measurement (A stands for the Arecibo Planetary 
Radar and G stands for the Goldstone Solar System Radar at DSS-14). 
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TABLE 3 

2000 ET70 HELIOCENTRIC ORBITAL ELEMENTS (SOLUTION 76) 



Element 



Value 



1-(T Uncertainty 



eccentricity 


0.123620379 


6.3 X 


10" 




semi-major axis (AU) 


0.9466347364 


1.2 X 


10" 


9 


inclination (degrees) 


22.3232174 


1.2 X 


10" 


6 


longitude of 










ascending node (degrees) 


331.16730395 


9.8 X 


10^ 


7 


argument of 










perihelion (degrees) 


46.106698 


1.1 X 


10" 


5 


Mean anomaly (degrees) 


84.37370818 


4.2 X 


10" 


7 



Note. — All orbital elements are specified at epoch 2012 
Dec 15.0 barycentric dynamical time (TDB) in the heliocen- 
tric ecliptic reference frame of J2000. The corresponding or- 
bital period is (336.41246710 ± 5.5 X 10"'') days. 



echo is given by: 



B 



cos S. 



(8) 



Here D is the diameter of the object, P is its apparent 
spin period, A is the radar wavelength, and S is the sub- 
radar latitude. As 6 increases, B decreases. In images 
obtained at similar rotational phases, the bandwidth ex- 
tent of the asteroid increased from February 12 to 20, 
indicating that our view was more equatorial towards 
the end of the observing campaign (Fig.jl]). For example 
in Fig. [8] images 10, 33, and 68 are at similar rotational 
phases and their bandwidths are ^--^3.7 Hz, ~4.9 Hz, and 
6.3 Hz, respectively. 

6. SPIN VECTOR 



We used the SHAPE software ( |Hudson||1993[ |Magri 
et al. |2007| to fit a shape model to the radar images 
and to estimate the spin vector of 2000 ET70. Since 
SHAPE is not particularly effective at fitting the spin 
axis orientation and spin period of the shape model, we 
carried out an extensive search for these parameters in an 
iterative manner. We performed two iterations each in 
our search for the spin axis orientation and spin period. 

Our initial estimate of the spin period came from the 
time interval between repeating rotational phases of the 
object captured in the images. Fig. [3] shows the object 
at similar orientations in images taken on different days. 
The time interval in two of those cases is ^72 hours and 
in the third case is ~45 hours. This indicates that the 
spin period of the object is close to a common factor of 
the two, that is, 9 hours or a sub-multiple of 9 hours. 
Fig. |4] shows two images taken 22.5 hours apart and the 
object is not close to similar orientations in these two 
images, ruling out all periods that are factors of 22.5 
hours. Thus we are left with a period close to 9 hours. 

Using an initially fixed spin period of 9 hours, we per- 
formed an extensive search for the spin axis orientation. 
The search consisted of fitting shape models to the im- 
ages under various assumptions for the spin axis orien- 
tation. We covered the entire celestial sphere with uni- 
form angular separations of 15° between neighboring trial 




Fig. 3. — 2000 ET70 range-Doppler images showing the asteroid 
at similar rotational phases on different days. The images were 
obtained at MJDs 55970.40 and 55973.38 (~72 hours apart) (top), 
55971.35 and 55974.33 (~72 hours apart) (center), 55969.46 and 
55971.34 (^45 hours apart) (bottom), suggesting a spin period of 
~9 hours or a sub-multiple of ~9 hours. In these images radar 
illumination is from the top, range increases towards the bottom, 
Doppler frequency increases to the right, and the asteroid spin 
results in counter-clockwise rotation. 

poles. At this stage we used a subset of the images in 
order to decrease the computational burden. We chose 
images showing sharp features that repeated on different 
days because they provide good constraints on the spin 
state. Images obtained on February 12 and 15 satis- 
fied these criteria and we used images with receive times 
spanning MJDs 55969.426 to 55969.463 and 55972.335 
to 55972.423 for this search. We started with triaxial el- 
lipsoid shapes and allowed the least-squares fitting pro- 
cedure to adjust all three ellipsoid dimensions in order 
to provide the best match between model and images. 
We then used shapes defined by a vertex model with 
500 vertices and 996 triangular facets. The fitting proce- 
dure was allowed to adjust the positions of the vertices 
to minimize the misfit. In all of these fits the spin axis 
orientation and the spin rate were held constant. Radar 
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Fig. 4.— 2000 ET70 images obtained at MJDs 55970.40 and 
55971.34, or 22.5 hours apart. The asteroid is not close to similar 
rotational phases in these images, ruling out spin periods that are 
sub-multiples of ~22.5 hours. A ^4. 5-hour spin period is therefore 
ruled out, leaving ~9 hours as the only plausible value. 

scattering parameters R and C described in Eq. ([t]) were 
allowed to float. We found the best shape model nt with 
the spin pole at echptic longitude (A) — 60° and ecliptic 
latitude (/?) = -60°. 

Using the best fit spin pole from the previous step, 
we proceeded to estimate the spin period with greater 
precision. We tried spin rates in increments of 2° /day 
from 960°/day (P = 9 hours) to 970°/day (P = 8.91 
hours) to fit vertex shape models with 500 vertices to 
the images. This time we used a more extensive dataset 
consisting of all images obtained from Arecibo. As in the 
previous step, only the radar scattering parameters and 
the shape parameters were allowed to float in addition 
to the parameter of interest. We found the best agree- 
ment between model and observations with a spin rate 
of 964°/day (Period = 8.963 hours). 

The second iteration of the spin-axis orientation search 
was similar to the first one except that we used a spin 
period of 8.963 hours and used the complete dataset con- 
sisting of all the Arecibo and Goldstone images from 
February. The Arecibo images from August were not 
used because of their low SNR and resolution. This pro- 
cedure was very effective in constraining the possible spin 
axis orientations to a small region of the celestial sphere 
(Fig. [5]) . We performed a higher resolution search within 
this region with spin poles ranging in A from 64° to 104° 
and (3 from 30° to 60° with step sizes of 4° in A and 
5° in p. For this step we performed triaxial ellipsoid 
fits followed by spherical harmonics model fits, adjusting 
spherical harmonic coefficients up to degree and order 
10. Our best estimate of the spin axis orientation is 
A = 80° and (3 = -50°, with 10° uncertainties. Shape 
models with spin axis orientations within this region have 
similar appearance upon visual inspection. We also at- 
tempted shape model fits with the best-fit prograde pole 
at A = 232° and (3 — 75° . We found that we were unable 
to match the observed bandwidths and ruled out the pro- 
grade solution. 2000 ET70 is a retrograde spinner, just 
like the majority of NEAs ( |La Spina et al. 2004). Our 
adopted spin pole (A — 80 , (3 = — 50°)^s at an an- 
gle of ^160° from the heliocentric orbit pole (A = 241°, 
P = 68°). 

We used the best spherical harmonics shape model 
from the previous search to perform a second iteration 
of the spin period search. This time we fit spherical har- 
monics shape models using spin periods in increments of 



0.001 hours from 8.940 hours to 8.980 hours. The reduced 
chi-squares of the shape models are shown in Fig. [6j We 
visually verified the quality of the fits and adopted a spin 
period of 8.960 ± 0.02 hours. The assigned uncertainty is 
conservative because a 0.02 hour difference in spin period 
amounts to an easily detectable ^25° offset in rotational 
phase over the primary observing window. 

7. SHAPE MODEL 

As a starting point for our final shape modeling ef- 
forts we used the results of the spin state determination 
(Section |6]), specifically the best-fit spherical harmonics 
shape model with a spin period of 8.960 hours and a 
spin pole at A = 80° and (3 = —50°. We proceeded to fit 
all the radar images obtained in February and OC CW 
spectra from Fig. |2] with a vertex model having 2000 ver- 
tices and 3996 facets. At this step the vertex locations 
and the radar scattering parameters were fit for, but the 
spin vector was held fixed. We minimized an objective 
function that consists of the sum of squares of residu- 
als between model and actual images, plus a number of 
weighted penalty functions designed to favor models with 
uniform density, principal axis rotation, and a reaso nably 
smooth surface ( Hudson|[T993{ |Magri et al.||2007[ ). The 
choice of weights in the penalty function is subjective, 
so the shape model solution is not unique. We tried to 
restrict the weights to the minimum value at which the 
penalty functions were effective. The minimization pro- 
cedure with our choice of penalty functions produced a 
detailed shape model for 2000 ET70 (Fig. |7] and Table [I]). 
The agreement between model and data is generally ex- 
cellent (Fig. |8] and Fig. [2]). The overall shape is roughly a 
triaxial ellipsoid with extents along the principal axes of 
^2.61 km, ^2.22 km, and ^2.04 km, which are roughly 
the same as the dynamically equivalent equal volume el- 
lipsoid (DEEVE) dimensions listed in table [4] 

The region around the north pole has two ridges that 
are 1-1.5 km in length and almost 100 m higher than 
their surroundings. These ridges enclose a concavity that 
seems more asymmetric than most impact craters. Along 
the negative x-axis a large protrusion is visible. Such 
a feature could arise if the asteroid were made up of 
multiple large components resting on each other. 

For this shape model, the best fit values for the radar 
scattering parameters, R and C, were 1.9 and 1.4, re- 
spectively (section |4]) . 

8. GRAVITATIONAL ENVIRONMENT 

We used our best-fit shape model and a uniform den- 
sity assumption of 2000 kg m~^, a reasonab le density for 
rubble-pile NEAs, (e.g., Ostro et "aL]|2006 ), to c ompute 
the gravity field at the surface of the asteroid (Werner 
and Scheeres 1997). The acceleration on the surface is 
the sum of the gravitational acceleration due to the as- 
teroid's mass and the centrifugal acceleration due to the 
asteroid's spin. An acceleration vector was computed 
at the center of each facet. Fig. [9] shows the variation 
of the magnitude of this acceleration over the surface of 
the asteroid. The acceleration on the surface varies be- 
tween 0.54 mm s~^ to 0.64 mm s~^, which is 4 orders of 
magnitude smaller than that experienced on Earth and 
2 orders of magnitude smaller than that on Vesta. Cen- 
trifugal acceleration makes a significant contribution to 
the total acceleration and varies from zero at the poles to 
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Fig. 5. — Contour plot showing goodness of fit (reduced chi-square) of shape model fits under various assumptions for the ecliptic longitude 
(A) and ecliptic latitude of the spin axis. We adopted a best fit spin axis orientation of A = 80° and /3 = —50° with a 10° uncertainty. 



TABLE 4 

2000 ET70 SHAPE MODEL PARAMETERS 



Parameters 




Value 


Extents along 


X 


2.614 ± 5% 


principal axes (km) 


y 


2.222 ± 5% 




z 


2.041 ± 5% 


Surface Area (km'^) 




16.68 ± 10% 


Volume (km'^) 




6.065 ± 15% 


Moment of inertia ratios 


AjC 


0.800 ± 10% 




B/C 


0.956 ± 10% 


Equivalent diameter (km) 




2.262 ± 5% 


DEEVE extents (km) 


X 


2.560 ± 5% 


y 


2.187 ± 5% 




z 


2.069 ± 5% 


Spin pole (A,/3) (°) 




(80, -50) ± 10 



Note. — The shape model consists of 2000 ver- 
tices and 3996 triangular facets, corresponding to an 
effective surface resolution of ^lOO m. The moment 
of inertia ratios were calculated assuming homoge- 
neous density. Here A, B, and C are the principal 
moments of inertia, such that A < B < C . Equiva- 
lent diameter is the diameter of a sphere having the 
same volume as that of the shape model. A dynam- 
ically equivalent equal volume ellipsoid (DEEVE) is 
an ellipsoid with uniform density having the same 
volume and moment of inertia ratios as the shape 
model. 
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u 0.386 
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Spin PefiocJ (hours) 

Fig. 6. — Goodness of fit (reduced chi-square) of degree-and- 
order 10 spherical harmonics shape models of 2000 ET70 as a func- 
tion of assumed spin period. We adopted a best fit spin period of 
8.960±0.02 hours. 



^-^0.049 mm s^^ (about 10% of the total acceleration) on 
the most protruding regions of the equator. For compar- 
ison, on Earth, centrifugal acceleration contributes less 
than 0.5% to the total acceleration at the equator. 

Fig. [To] shows the gravitational slope variation over 
the asteroid's surface. The gravitational slope is the an- 
gle that the local gravitational acceleration vector makes 
with the inward pointing surface-normal vector. The av- 
erage slope is 9.5° and it reaches up to ^^33° on the sides 
of the ridges at the north pole. The sides of the ridges 
facing the north pole are steeper than the opposite sides. 
Less than 1% of the surface has slopes greater than 30°, 
the approximate angle of repose of sand, indicating a 



2 km 




view from - x 



view from -y 



view from -z 



Fig. 7. — Best-fit vertex shape model of 2000 ET70 as seen along the three principal axes x, y, and z. For principal axis rotation the spin 
axis is aligned with the z axis. Yellow regions at the south pole have radar incidence angles > 60° and hence are not well constrained. The 
shape model has 2000 vertices and 3996 triangular facets giving an effective surface resolution of ~100 m. 



relaxed surface. The regions near the north pole where 
the slopes are the highest may be composed of competent 
rocks exposed at the surface. Along-the-surface accelera- 
tions on the sides of the ridges are as high as 0.34 mm s~^ , 
so one might expect to see accumulated debris in the val- 
ley between the two ridges. 

9. DISCUSSION 



rain is playing an important role, there are albedo vari- 
ations over the surface of the asteroid, or a combination 
of these factors. 

also 



Alvarez et al. 



Whiteley| ( |2001[ ) estimated 2000 ET70's absolute mag- 
nitude [H) to be 18.2. The geometric albedo (pv) of an 
asteroid i s related to its effective d iameter (D) and its H 
value by ( Pravec and Harris||2007 ): 

^-0.2H-\ 2 



Pv 



1329 km X 10" 



D 



(9) 



Using |Whiteley| ( |2001 )'s ir=18.2, the above equation 
yields py=0.0iy±0.002 for an asteroid with an effective 
diameter of 2.2 km ± 6% (Fig. 11). This uncertainty 
on Pv is due to the diameter unc ertainty only as n o un- 
certainties on H were quoted in [Whiteleyl ((2001|. Al- 
though py = 0.019 is possible, it seems unreasona bly low 
compared to the albedos of other NBAs (Thomas et al.^ 
2011 Stuart and Binzel||2004 ). More conimon values of 
Pv would require lower values of the absolute magnitude 
(e.g., ^^=0.04 requires H—17.5). The range of possible 
H and pv values that are con sistent with the radar size 
est imates is shown in Fi g. |11| 

Alvarez et al. ( 2012[ ) observed the asteroid between 



February 19 to 24, when the view was close to the equa- 
tor, and reported a lightcurve amplitude of 0.60 ± 0.07 
mag. If the asteroid was approximated by a triaxial ellip- 
soid with uniform albedo, this amplitude would suggest 
an elongation (ratio of equatorial axes) approximately 
between 1.28 and 1.36. Our shape model indicates that 
this ratio is ^1.18, suggesting that either the ellipsoid ap- 
proximation is poor, the lightcurve amplitude is on the 
lower end of the range above, shadowing due to the ter- 



(|2012|) also report a period of 
8.947 ± 0.001 hours, which they refer to as the "syn- 
odic period", and which we refer to as the "lightcurve 
period" . The latter is a function of the intrinsic spin 
state of the asteroid and of the relative motion between 
the asteroid, the observer, and the Sun. Our shape mod- 
eling process relies on the intrinsic or sidereal period of 
the asteroid, which we estimate to be 8.960 ± 0.02 hours. 
This value can be used with our best-fit spin axis orien- 
tation and ephemerides to evaluate what we refer to as 
the "synodic period", i.e., the period that combines the 
(fixed) intrinsic rotation and the (variable) apparent ro- 
tation due to sky motion, irrespective of the position of 
the Sun. On Feb 19.0, the nominal synodic period was 
8.937 hours, whereas on Feb 25.0, the nominal synodic 
period was 8.943 hours. These synodic periods are close 
to the reported lightcurve period, but cannot be directly 
compared to it as they measure slightly different phe- 
nomena. 

Arecibo and Goldstone radar observations of 
2000 ET70 allowed us to provide a detailed characteri- 
zation of a potentially hazardous asteroid, including its 
size, shape, spin state, scattering properties, and grav- 
itational environment. These techniques are applicable 
to a substantial fraction of known NEAs that make 
close approaches to Earth within ^0.1 AU. Radar-based 
physical properties for this and other asteroids are 
available at http://radarastronomy.org 
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Fig. 8. — Comparison between radar rangc-Doppler images, corresponding synthetic radar images generated using our best-fit shape 
model, and plane-of-sky (POS) projections of the shape model. Range-Doppler images are oriented as in Fig.ls] and the POS projections 
are oriented north-up. Time increases from top to bottom within each panel and from left to right. The spatiaTresolution of the images is 
given in Table ^ 
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Fig. 9. — This figure shows the magnitude of the vector sum of acceleration due to gravity and centrifugal acceleration computed at the 
center of each facet of the shape model. We assumed a uniform density of 2000 kg m""^ and used our measured spin period value of 8.96 
hours. 
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Fig. 10. — This figure sliows ttie gravitational slopes computed at the center of each facet of the shape model. We assumed a uniform 
density of 2000 kg m~^ and used our measured spin period value of 8.96 hours. 
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Fig. 11. — Absolute magnitude vs. geometric albedo for an effec- 
tive diameter of 2.2 km (dark line) and a diameter uncertainty of 
6% (light lines). The horizontal dashed l ine corresponds t o the ab- 
solute mag nitude (H=18.2) reported by'Whiteley| | |2001| l, and the 
vertical dashed lines indicate the corresponding geometric albedos. 
These albedo values are unusually low, suggesting an extraordinar- 
ily dark object or an error in the reported H value. 
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